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Abstract. 

The field of warp research has been dominated by analytical methods to investigate 
potential solutions. However, these approaches often favor simple metric forms that 
facilitate analysis but ultimately limit the range of exploration of novel solutions. So far 
the proposed solutions have been unphysical, requiring energy condition violations and 
large energy requirements. To overcome the analytical limitations in warp research, 
we introduce Warp Factory: a numerical toolkit designed for modeling warp drive 
spacetimes. By leveraging numerical analysis, Warp Factory enables the examination of 
general warp drive geometries by evaluating the Einstein field equations and computing 
energy conditions. Furthermore, this comprehensive toolkit provides the determination 
of metric scalars and insightful visualizations in both 2D and 3D, offering a deeper 
understanding of metrics and their corresponding stress-energy tensors. The paper 
delves into the methodology employed by Warp Factory in evaluating the physicality 
of warp drive spacetimes and highlights its application in assessing commonly modeled 
warp drive metrics. By leveraging the capabilities of Warp Factory, we aim to further 
warp drive research and hopefully bring us closer to realizing physically achievable 
warp drives. 
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1. Introduction 


Warp drive research has been mostly driven by analytic studies of various metrics. Such 
limitations have led to a narrow perspective when it comes to warp drive physicality, 
especially due to the non-linear nature of the gravitational field equations. We developed 
Warp Factory to broaden the view of warp drive research. Our numerical toolkit can 
evaluate the Einstein equations to facilitate the analysis of the stress-energy tensor and 
its consequent energy condition satisfaction or violation. 

The first warp drive solutions developed by Alcubierre and then further 
developed by Van Den Broeck [16] enabled passengers inside the warp bubble to move 
at arbitrary speeds but at the cost of both enormous energy density and unphysical 
stress-energy tensors using negative energy. The recent solutions from Lentz [II] and 
Fell [5| both find solutions with positive Eulerian energy density. However, it has been 
argued that these still both violate the weak energy condition, as shown generally by 
and by Fell himself. Since the publication of [2| an extended class of possible solutions 
has been proposed for warp drives, with an emphasis on physical solutions which require 
subluminal speeds. 

The key challenge in warp research is two-fold: The first is evaluating the 
complicated Einstein field equations to find stress-energy tensors and the second is 
in evaluating those stress-energy tensors for physicality. Proving physicality is a unique 
focus in the field of warp drive research that is often not required when solving general 
astrophysical problems. For simple warp drives, such as those with symmetry and 
time-invariance, analytical solutions for the stress-energy tensor can usually be easily 
determined. Once the analytic stress-energy tensor is known, the physicality of the 
solution is then evaluated by contracting the stress-energy tensor with the velocity 
of different observers. Often the contraction process to show physicality is incorrectly 
limited to a single, easily-defined observer. For simple metrics, this is typically sufficient 
to prove that a solution is not physical. In a general case, physicality is more difficult 
to prove and single observers are not sufficient. The broader cases of warp solutions, 
such as those with Schwarzschild boundaries, asymmetries, and accelerations, will have 
complicated metrics that are inaccessible to analytical approaches. This presents a 
significant challenge to the development of new warp solutions. 

To overcome the difficulties of creating and analyzing general warp drives, Warp 
Factory) a numerical analysis toolkit, has been developed that provides a general 
numerical evaluation of the Einstein field equations and physicality of the stress-energy 
tensor across a range of observers for arbitrary input metrics. Along with testing for 
physicality, the toolkit can provide other analyses, such as constructing metric scalars 
and analyzing momentum flows. This numerical approach offers robust methods to 
evaluate physicality and provides the ability to visualize the metric and stress-energy 
tensors in order to gain valuable insights. 

In this paper, we establish the methods used in Warp Factory to construct and 


5 https://github.com/NerdsWithAttitudes/WarpFactory 
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analyze spacetimes and demonstrate this analysis on several example warp drive metrics. 
Warp theory and the definition of physicality are introduced in Section [2] The evaluation 
of the Einstein field equations, the methods used in conducting stress-energy tensor 
analysis, and the solution of metric scalars are explained in Section [8| Then, in Section 
these methods are applied to evaluate metrics commonly considered in warp literature. 
Finally, Section [5] summarizes and discusses the results, and Section [6] concludes with 
closing remarks. 


1.1. Key Results 


(i) Numerical methods are established and demonstrated for evaluating the physicality 
of metrics using a custom framework in MATLAB called Warp Factory. 


(ii) The Alcubierre, Van Den Broeck, Bobrick-Martire Modified Time, and Lentz- 
inspired warp solutions are evaluated, showcasing their stress-energy components, 
energy condition violations, and metric scalars. 


2. Warp Theory 


The field of warp mechanics does not have a unified definition of what is required by 
a warp solution. However, there are several aspects that emerge as important features 
of any solution. A warp drive aims to provide a controlled and comfortable journey 
through space. This goal motivates the three core features of a warp drive: geodesic 
transport, a flat passenger region, and a spatially compact and comoving warp bubble. 


(i) Geodesic transport. Warp drives should transport passengers from point A to point 
B along a geodesic trajectory. This means the passengers inside the warp drive 
do not experience local acceleration while being transported"| For a non-trivial 
solution, the passengers should not ‘already be going to’ point B. One possible 
scenario is that passengers start at rest (relative to points A and B) at point A, 
are transported to point B, and are then, again, at rest relative to both points. 
Additionally, this transport should occur for a limited region in spacetime. 


(ii) Empty passenger region. There should exist a defined passenger volume that is 
both vacuunt’| (T“” = 0) and free from large tidal forces. 


(iii) A spatially bounded, comoving bubble. The stress-energy required for the geodesic 
transport should not extend to infinityf] The stress-energy of the bubble also needs 
to move along with the transported observers. The requirement of it comoving with 
passengers distinguishes warp drive solutions from Krasnikov tubes [4]. 


6 In the case that a local acceleration is desired, for example, 1g, the statement becomes the local 
acceleration of passengers should be limited. 

T Ignoring the mass of the passengers. 

8 Perhaps some energy could be radiated to infinity but that energy should be causally connected to 
the bubble. 
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This framework establishes the primary goals of any generic warp solution. The 
warp metrics proposed in the literature have all satisfied different features in this list. 
For example, Alcubierre has a flat passenger volume and a warp bubble whose spatial 
extent is determined by a shape parameter) 


2.1. Conception of Warp 


In Star Trek and other science fiction, the notion of a warp drive has several distinct 
features which may or may not apply to real, physical, warp drives including faster- 
than-light movement, self-acceleration, a normal rate of time-passage, and possessing a 
bubble of energy that surrounds a ship. 


e Moving Faster-than-light: The creation of warp has historically been motivated by 
the ability to potentially travel at superluminal speeds, but the broader solutions of 
warp drives are not restricted to this criteria and interesting solutions for subluminal 
speeds also exist. In the recent paradigm by Bobrick [2], warp drives are classified 
into regimes of solutions based on speed. Likely, physical solutions are restricted 
to the subluminal regime [2]. 


e Self-acceleration: To transport observers at rest to some other point in spacetime 
requires a coordinate acceleration as measured by an external observer. For objects 
in flat space, acceleration requires the transfer of momentum of the accelerating 
system. Within a curved spacetime, observers can travel along geodesics which 
can be perceived by other observers as a relative acceleration. In addition, the 
constraint of a comoving bubble also requires the bubble itself to self-accelerate or 
be conventionally propelled. 


e Passenger Time: In line with the Star Trek conception of warp drives, the first 
warp solution from Alcubierre was constructed such that the passengers and 
external observers both experience the same passage of time (proper time). This 
construction exists in most of the solutions developed. However, allowing the 
modification of the lapse rate might be an important part of physical solutions. 
Warp shells with a bubble consisting of regular matter, such as a matter shell, may 
require a change in lapse to the passenger volume [7]. 


e Spaceships and Warp Bubbles: In Star Trek, the ship creates some kind of radiated 
energy field that dissipates when the warp effect is removed. This does not need 
to be the case in general. For example, a bubble might be constructed by a shell 
of matter whose stress-energy tensor changes during the warp journey but returns 
back to a regular matter shell when the warp effect is removed. 


9 While the functions used by Alcubierre are not exactly zero at r >> Rbubble, it is, for practical 
purposes, approximately compact. This is also true for the passenger volume flatness. 
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2.2. Defining Physicality 


The physicality of warp metrics is established through a set of conditions, known as 
energy conditions, which are imposed on the stress-energy tensor components [10]. 
Among these conditions, the most frequently used and implemented within Warp 
Factory are: 


(i) The Point-wise Null Energy Condition (NEC) is the null (light-ray) frame’s energy 
density observation, at all points in spacetime. For this condition to be physical, 
the density should be non-negative for any null observer at every point in spacetime. 


(ii) The Point-wise Weak Energy Condition (WEC) is a timelike frame’s energy density 
observation, at all points in spacetime. For this condition to be physical, the density 
should be non-negative for any timelike observer at every point in spacetime. 


(iii) The Point-wise Strong Energy Condition (SEC) is a measure of the experienced 
tidal effect from the matter present acting on timelike observers. For this condition 
to be physical, tidal effects must be non-negative, meaning that matter should 
gravitate together at every point in spacetime. 


(iv) The Point-wise Dominant Energy Condition (DEC) is an observation of the matter 
flow rate at all points in spacetime. For this condition to be physical, the velocity 
must be observed to be less than the speed of light. 


All known forms of classical matter obey these conditions. However, quantum 
systems can violate them for local points [6]. Therefore, a softer version of the NEC, 
WEC, and SEC is the average of each of those conditions where each energy condition is 
integrated over a geodesic path instead of evaluated pointwise. In addition, there is an 
argument to be made to drop the SEC since dark energy and inflation appear to violate 
it [9]. For this analysis, however, we will evaluate the more stringent requirements of 
satisfaction for all energy conditions for all observers, everywhere. 

Satisfaction of the energy conditions alone cannot guarantee the prospect of 
manmade warp drives. Any warp engineer will have to include practical requirements 
for the drive, such as: 


e Conservation of stress-energy V T”” = 0. This is the conservation law of general 
relativity. If the stress-energy tensor is given by (11), then, by construction, any 
metric that is continuous up to the second order will satisfy this condition per the 
contracted Bianchi identities [8]. 

e No horizons of inaccessibility between the passenger volume and the external 
spacetime. Warp structures should always be limited to a regime where observers 
can be connected to each other across the warp bubble and to the external space. 

e No violation of causality. Closed time-like curves are undesirable. 

e Reasonable total energy and energy densities Any practical warp bubble should use 


an attainable amount of matter in configurations of densities that are practically 
achievable. 
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2.8. Observing Stress-Energy 


The requirements of a warp drive are a description of the spacetime curvature given 
by the metric tensor guu. However, to find the energy, momenta, pressures, and shears 
required to generate that curvature in spacetime, the Einstein field equations must be 
used to compute the stress-energy tensor from the metric tensor. Einstein notion will 
be used for summations where the 0 component of an indexed object will be reserved 
for a time-like direction whereas the 1-3 components are for space-like directions. Latin 
indices denote summations from one to three and Greek indices denote summations 
from zero to three. 

The components of the stress-energy tensor represent the coordinate observer’!)] 
perspective on the stress-energy throughout the space. To better understand the actual 
distribution of energy, momenta, pressures, and shears, local observations of these 
quantities must be made using observers at those same points in spacetime. 


2.8.1. Constructing Observers Observers in general relativity are fully specified by a 
timelike vector and three spatial vectors, with the timelike vector usually specified as a 
four-velocity. This can also be put into a tetrad formalisn{™] where each observer has a 
tetrad that is associated with their frame. The tetrad can be given as a matrix whose 
columns are constructed from the time vector or four-velocity e5 of the observer and its 
spatial vectors e: 


eË e5 ek (1) 


When evaluating the field equation from the metric tensor what is returned is the stress- 
energy as represented in the coordinate frame. Generally, we can express any frame’s 
stress-energy tensor Tap in terms of the coordinate tensor and the frame’s tetrads: 


Tho = Twejeo (2) 
where the hats represent the observer frame. If the observer is aligned with the frame 


defined by the tetrad, then the observables of energy (p), momenta (p), pressures (P), 


and shears (g) are directly given as: 
p= Tron!n® == 00 
Pi =— non? = —Íġğ = he (3 
P; = Toti ti = Ta =T” 


a = Ty a 


10 The coordinate observers represent observers whose local metric is Minkowskian in the specified 
metric. For example, for the Schwarzschild metric in Schwarzschild coordinates, the coordinate 
observers are stationary observers at infinity. 

11 Often observers are only associated with a four-velocity, however, this leaves the spatial vectors 
undefined. Fully defined observers with a spatial vector will have a unique tetrad. For many problems, 
only the four-velocity is required. 
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where, i # j and, in this specific instance, n” and y# are defined as aligned with the 
coordinate system, resulting in: 


0 0 0 0 

h â 0 1 0 0 
# = (1,0,0,0 = 4 
n (1,0,0,0), Vi 001 0 (4) 

0 00 1 


2.3.2. Eulerian Observers There is a definition of a particular tetrad that can be 
expressed in terms of the metric tensor: 


g” = ehen”? (5) 


With this tetrad, we can construct the inverse of this statement where the tetrad acts 
on the covariant form of the metric to return the locally-Minkowskian metric: 


Juv nes = Naw (6) 


Finding the observer tetrad is then simply solving this equation given the metric tensor. 
This solution is not unique unless the temporal and spatial orientation is first defined, 
which depends on the selected observer. 

A special set of observers is the set of Eulerian Observers which is defined as having 
a tetrad with a four-velocity that is orthogonal to the spatial hypersurface (not moving 
in space over time). Using 3+1 formalism the foliation of the spacetime is constructed 


as [8]: 
-a° + BB Bj 
= T 
ad í pi 7 ( ) 


where a represents the lapse rate and where the spatial metric Ņ;; raises the shift vector 


index ĝ; as: 
BY =” Bb; (8) 
We can then decompose any metric defined in the standard form into its 3+1 
components: 
Bi = Joi 
Vig = Jij (9) 


Q&Q = vy g” goi — Joo 


This formalism allows for a convenient definition of the four-velocity which is orthogonal 
to the spatial hypersurface: 


1 

eh =n" = & (1, —', — 6, — 6°) > Rhy = (=a@,0,0,0) (10) 
and the spatial terms of es can be selected freely assuming they are orthogonal to the 
time vector and each other. In Warp Factory we define the spatial vectors by selecting 
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that ĉ is orthogonal to the coordinate y and z surface and the remaining spatial vectors 
being orthogonal to each other (full solution and vector construction are detailed in 


Appendix AÌ**| With this setup, we have fully defined an Eulerian tetrad with a selected 


spatial orientation which can be solved symbolically from only the covariant metric. 


3. Methods of Warp Factory 


In Warp Factory, each of the 10 independent metric components (gv) can be specified 
at each point on the 4D spacetime grid. For correct evaluations, the metric is required 
to have a Minkowskian signature of (-,+,+,+) at all points. 


3.1. Evaluating the Field Equation 


The stress-energy tensor is determined from the metric through the field equations: 


ER ME (11) 
Aoa one 


In Warp Factory, the Ricci tensor R,, is calculated directly in terms of the metric 
derivatives: 


oe Pow _ Ogas _ gus P's \ 9p 
Rw = 2 2 (e T Orð” Ox ôr Ortðre) F 


aœ, ß=0 
3 
1 1 Jay Oggi ONA O9v5 Onan Ogue B y8 
5 - i 12 
t3 2 c Aa" Ôx” ' Ore Ox® xe Orë)’ I (12) 
a,8,7,d=0 
3 
1 Oguy Juy Jw gg Jag 5 
== _ 2 — aB ay 
M 2A ( ðr” Or” Ox dre Og J FF 


The Ricci scalar R is: 
sig (13) 


The metric derivatives are found using a fourth-order central finite difference method 
[I5] with a 1-meter grid spacing. The boundary points of the metric are cropped in the 
evaluations. See appendix [Appendix B] for details on accuracy. 

For the computation of T””, tetrads are selected to fulfill the requirements of (LO). 
A full solution in terms of the covariant metric terms are detailed in 
This alternative approach is used within Warp Factory as the metrics are not strictly 
constructed in 3+1 formalisn{"| 
12 This choice of alignment is arbitrary but needed to provide a unique solution. 


13 The metrics created and evaluated with this method must be able to be constructed in 3+1. The 
tetrad creation is equivalent 
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3.2. Evaluating Physicality 


In Warp Factory, to establish if a solution is physical, we evaluate the standard energy 
conditions (NEC, WEC, SEC, and DEC) for all points in discretized spacetime. All 
evaluations of the energy conditions require a contraction of the stress-energy tensor 
with an observer vector field that is either null or timelike. Numerically, this process 
involves sampling observer velocities and contracting them with the stress-energy tensor 
at all spacetime points on the coordinate grid. The stress-energy tensor and the energy 
conditions will be evaluated in the tetrad frame, T“”, which, at all points, is associated 
with a local Minkowskian metric. Each observer is defined by a four-velocity that can 
be either null or time-like. The evaluation of energy conditions does not require a 
specification of the spatial vectors of the tetrad if the tensor is transformed into any 
Eulerian frame. It can simply be represented as the contraction of the stress-energy 
tensor with four-velocities, so we drop the strict tetrad description here for convenience. 


3.2.1. Observer Vectors and Vector Fields The null vector for a given observer is kÊ, 
which, in a Cartesian locally-Minkowskian space, is defined as: 


kê = (1,@) (14) 


where w is some vector representing the spatial velocity of an observer that satisfies the 
condition of |v] = 1. For the set of null observers On, the direction of the spatial velocity 
is sampled from a set of vectors that map to evenly distributed points on a sphere, which, 
in a spherical coordinate representation, result in the set of four-velocities that map to 
each observer as: 


O! + kê = (1, sin(6,) cos(¢,), sin(@;) sin(¢,), cos(@,)) 
On = (15) 
O” — kê = (1, sin(@,) cos(¢,), sin(@,) sin(dn), cos(On)) 


The resulting vectors for each observer O” are then normalized. For timelike 
vectors, VÊ, the four-velocity is defined in the same manner but |t| is additionally 
scaled from 0 to 1 for a specified number of timelike samples for each of the spatial 
velocity direction samples. The timelike observers, O, are constructed as: 


O! + V° = (1, sin(,) cos(¢1)s1, sin(@;) sin(¢1)s1, cos(41)s1) 
ae (16) 
O” + VÊ = (1, sin(On) Cos(bn)Sn; SiN(An) Sin(Gn) Sn, COS(An) Sn) 


where 0 < sp < 1. All vectors are then normalized. 

Since the metric, gap, is locally-Minkowskian everywhere, the same vector can be 
contracted at all spacetime points X for Tao(X). This means that kê and V° can be 
treated as uniform vector fields since their values do not vary over the spacetime. 
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3.2.2. Energy Conditions The Null Energy Condition (NEC) is given by the 
contraction of stress-energy tensor with all observer null vector fields: 

En(X) = Tho(X)kAk? > 0 V k’ (17) 


Similarly, the Weak Energy Condition (WEC) is a contraction of the stress-energy tensor 
with all observer timelike vector fields: 


Ew(X) = Tel X) VV? > 0 V VA (18) 
The Strong Energy Condition (SEC) is found by the equation: 
1 "ORE: 
E5(X) = (Ta) = ITD) VV” > OV yr (19) 


Lastly, the Dominant Energy Condition (DEC) is satisfied when both the WEC is 
greater than zero and when the contracted mixed tensor with a timelike vector field is 
less than zero. The latter part of the DEC is given by: 


T#(X) = —T#,(X)V? (20) 


where this new vector field TÊ (X) must be future pointing and either timelike or null. 
This is found by the equation: 


Ep(X) = nao T*(X)T?(X) <0 V VA (21) 


For easy comparison to other energy conditions, the square root of this value is taken 
while flipping the overall sign: 


Ep(X) = —sign(Ep(X)) V'lEp(X)| (22) 

Therefore, just like the other energy conditions, when the DEC is evaluated, a negative 

result represents the violation of the condition. In Warp Factory, only the evaluation of 
=Ep(X) is returned as the dominant energy condition. 

For each energy condition, the minimum value across all observers is taken for each 


point in the spacetime grid. These minimum values are what is plotted in the energy 
condition evaluations presented in Section [4] 


3.8. Momentum Flowlines 


Warp metrics generate terms beyond just the energy density. One unique feature that 
often exists in warp drive spacetimes is momentum flux p;, which describes the motion 
of the energy density. 

This flow can be visualized by constructing a set of flowlines that are built from 
the vector field of the momentum density, given by Q: 


U(X) = (P(X), po(X), ps(X)) (23) 


Each of the flowlines is found by propagating a virtual particle from a starting point 
through the vector field for a given number of steps. It should be noted that while this 
provides insight into the structure of the momentum density, these lines are not those 
of matter moving along geodesic paths. 
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3.4. Metric Scalar Analysis 


The metric tensor can also be understood by inspecting its associated geometric scalars, 
which are constructed from the metric tensor by projections onto a timelike vector field 
U". Any timelike vector field can be used, but the specific vector field assumed in Warp 
Factory is that of a four-velocity orthogonal to the spatial hypersurface as given by 
Ur =n. 


3.4.1. Expansion Scalar The expansion scalar describes how the volume of a region 
changes. We define this scalar 0 as: 


0 = g 0w 24 
H 


where ĝu» is the stress tensor, which is computed using the projection tensor P% that 
operates on the covariant derivative of the vector field U# aq'4] 


Ou = P PIV (Ua) (25) 
and the projection tensor is given by: 


Pw = Guv + UU, (26) 


3.4.2. Shear Scalar The shear scalar describes the relative stretching while preserving 
the object’s volume. The shear scalar g? is defined as: 


1 
= ae (27) 


where „v is the shear tensor, which is similarly computed using the stress tensor along 
with the projection tensor P,,,, shear tensor 6,,,, and shear scalar 0 as: 


= yp ee ws (28) 


3.4.8. Vorticity Scalar The vorticity scalar is the magnitude of the twist in spacetime, 
given by: 


1 
w? = gouw” (29) 
Where wp is the vorticity tensor, which is computed using the projection tensor as 
wig =P P Vaa (30) 


Given a vector field that is orthogonal to the spatial hypersurface, the vorticity scalar 
will always be zero, but this is not generally true for an arbitrary vector field. 


14 Note: X (ab) = Xab + Xba and X [ab] = Xab — Xba- 
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4. Metric Evaluations 


The numerical approach that Warp Factory enables provides greater insight into the 
details of warp spacetimes. As a demonstration, we will consider several popular metrics 
discussed in the literature. 


4.1. Alcubierre Metric 


Alcubierre introduced the first general relativistic warp drive in 1994 [i]. By adding 
a localized and spherical shift to a Minkowski spacetime, one obtains a flat passenger 
volume moving along a pre-defined curve. The Alcubierre metric is given as: 


ds? = —dt? + (dx — v,f (rs) dt)? + dy? + dz? (31) 


where 


_ tanh (o (rs + R)) — tanh (ø (rs — R)) 
2 tanh(o R) (32) 


rs = Ve- 20 — vst)? + (y — yo)” + (z — 20)" 


f (rs) 


To analyze the metric, the commonly used Galilean transformation is applied: 


x = T — Vst 


t=t se) 


This makes the metric time-independent and amounts to an addition of the shift vector: 


B, = Bo + Us (34) 


The metric being in a comoving frame removes the necessity for time derivatives and 
reduces the computation required. The Alcubierre metric violates all energy conditions 
within the bubble region, as it necessitates negative energy density, as Alcubierre 
demonstrated. This fact is confirmed through the analysis conducted using Warp 
Factory. The metric solved used the following parameters: v, = 0.1 c, R = 300 m, 
o = 0.015 m™t, x = yo = % = 503 m. The ratio of R to sigma was chosen so 
that there is a substantial region of flat space inside of the passenger area and so that 
the derivatives of the shift were not extremely large across the bubble boundary. In 
addition, R was chosen sufficiently large to have high accuracy when evaluating the 
stress-energy tensor and the energy conditions at a l-meter grid spacing. This creates a 
metric which is shown in Figure[I| The resulting stress-energy tensor is shown in Figure 
which depicts the energy density, and Figure |3| which presents the complete stress- 
energy components. Finally, Figure |4| showcases the energy conditions, and Figure 
shows the metric scalars. 
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Figure 1. Alcubierre metric comoving shift vector for the X-direction. The motion 
of the bubble is in the +X direction. The slice is for z = z, vs = 0.1 c, R = 300 m, 
o = 0.015 m7!, £o = yo = zo = 503 m. 
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Figure 2. Energy density for the Alcubierre metric. The motion of the bubble is in 
the +X direction. The slice is for z = zo, vs = 0.1 c, R = 300 m, o = 0.015 m7}, 
Lo = Yo = Zo = 503 m. Units are in J/m?. 
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Figure 3. The stress-energy tensor for the Alcubierre metric. The 700 component is 
shown in Figure [2] The motion of the bubble is in the +X direction. The slice is for 
z = 2%, vs = 0.1 c, R = 300 m, ø = 0.015 m™!, £o = yo = zo = 503 m. Units are in 
J/m?. 
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Figure 4. Energy condition evaluation for the Alcubierre metric. The motion of the 
bubble is in the +X. The slice is for z = zo, vs = 0.1 c, R = 300 m, o = 0.015 m7, 
zo = Yo = Zo = 503 m, Nobservers = 1000. The minimum value among all observers is 
shown. Units are in J/m?. 
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Figure 5. Expansion and shear scalars for the Alcubierre metric. The motion of the 
bubble is in the +X. The slice is for z = zo, vs = 0.1 c, R = 300 m, ø = 0.015 m™!, 
Lo = Yo = Zo = 503 m. 
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4.2. Van Den Broeck Metric 


Shortly after Alcubierre’s original paper, Van Den Broeck formulated his own variation 
of the Alcubierre metric. His approach involved employing two concentric regions of 
varying shift-vector and expansion terms. The objective was to diminish the need 
for negative energy through the increase in the volume of the passenger region. This 
formulation, as defined by Van Den Broeck’s initial publication, is given by: 


ds” = —dt? + B? (r,) [(da =v] (rs) dt)” + dy? + d2?| (35) 
where f(r.) is: 


f(r) =1for r< R 
O0<f(r,)<lfor R<r,<R+A (36) 
f(r, =0 for AAAs, 


and B(r,) is: 


B(r,)=l+afor r,<R 
1<B(r,.)<l+afor R<r,<R+A (37) 
Bre) = 1 for R+A<r, 


Here R is the inner region radius, R is the outer radius of the warp bubble, and A 
and A are the thickness of the transition regions. As in the Alcubierre metric analysis, 
the Galilean transformation to a constant velocity comoving frame will be used. The 
parameters used are v, = 0.1 c, a = 0.5, R = 350 m, R= 200 m, A= A= 40 m, 
To = Yo = Zo = 503 m. The metric component cross-sections are shown in Figure [6] 

Figure |7| illustrates the energy density profile, which shares a resemblance to 
Alcubierre but features a positive region within the transition into the passenger volume 
as it expands. Similarly, Figure p] demonstrates that this metric also violates energy 
conditions, similar to Alcubierre, but with rings and a non-violating region where the 
positive energy exists. The presence of negative energy in the Eulerian frame indicates 
the existence of a violation, but understanding the distinct rings of violation requires 
consideration of the remaining stress-energy tensor components, depicted in Figure 
The interplay between pressure and momentum density, particularly in areas of low 
energy density, gives rise to the primary regions of violation. Lastly, the Van Den 
Broeck metric scalars, shown in Figure have many similar features to Alcubierre, 
just divided along its concentric shell regions. 

Van Den Broeck’s approach does not address the key issue which creates the energy 
condition violations in the Alcubierre metric, but it may provide a unique method to 
reduce the total violation amount per passenger volume, thus creating a more efficient 
warp drive through its use of volume expansion. 
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Figure 6. Metric components for the comoving Van Den Broeck metric. The Van 
Den Broeck metric uses a shift vector in the X direction and expansion in 7;; terms. 
The motion of the bubble is in the +X direction. The slice is for z = zg, vs = 0.1 c, 
a = 0.5, R = 350 m, R= 200 m, A = Ã = 40 m, £o = yo = Zo = 503 m. 


Analyzing Warp Drive Spacetimes with Warp Factory 19 


roo 
1000 1.5 
900 i 
800 
0.5 
700 
470 
600 
= 500 J3 
> 
400 -1 
300 
-1.5 
200 
-2 
100 
-2.5 
0 


O 100 200 300 400 500 600 700 800 900 1000 x10% 
X [m] 


Figure 7. Energy density for the Van Den Broeck metric. The motion of the bubble 


is in the +X direction. The slice is for z = zo, vs = 0.1 c, a= 0.5, R = 350 m, R = 200 


m, A = A = 40 m, 20 = yo = zo = 503 m. Units are in J/m?. 
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Figure 8. Stress-energy tensor for the Van Den Broeck metric. The 700 component 
is shown in Figure The motion of the bubble is in the +X direction. The slice 
is for z = zo, vs = 0.1 c, a = 0.5, R = 350 m, R = 200 m, A= A= 40 m, 
£o = Yo = zo = 503 m. Units are in J/m?. 
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Figure 9. Energy conditions for the Van Den Broeck metric. The motion of the 


bubble is in the +X direction. The slice is for z = zp), Vs = 0.1 c, a= 0.5, R = 350 m, 


R = 200 m, A= A= 40 m, zo = yo = Zo = 503 m, Nobservers = 1000. The minimum 
value among all observers is shown which can be positive. Units are in J/m°. 
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Figure 10. The expansion and shear scalars for the Van Den Broeck metric. The 
motion of the bubble is in the +X direction. The slice is for z = zg, vs = 0.1 c, a = 0.5, 
R= 350 m, R= 200 m, A= A= 40 m, zo = yo = 2% = 503 m. 
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4.8. Bobrick-Martire Modified Time Metric 


Bobrick and Martire explore various alterations to the Alcubierre metric in their 
publication titled “Introducing Physical Warp Drives” [2]. Among these modifications, 
they specifically examine the incorporation of a dynamic lapse rate, which is detailed in 
section 4.5 of their work. The line element, where p and @ are cylindrical coordinates 
about the axis of direction of travel x, are given as: 


ds? = — ((1— f(r.))dt + A(r.)“1f dt)” + (dz — f (rs)us dt)? + dp? +p? d6? (38) 
where A(r,) is: 
ae K (r)>1 for r< R (39) 
1 fr r> R 


The inclusion of A(r,) changes the lapse rate inside the drive, hence the name “modified 
time metric”. Again, this metric is analyzed in the constant velocity comoving frame. 
In this evaluation, A(r,) is constructed with the same fall-off function and ø as the 
f(r.) function in Alcubierre. The parameters for the metric are v, = 0.1 c, Amaz = 2, 
R = 300 m, o = 0.015 m™!, £o = yo = % = 503 m. The metric components for this 
warp bubble are shown in Figure 


4000 Lapse Rate 4000 X Shift 
800 800 
— 600 8 = 600 
E E, 
> 400 7 >% 40 
200 0.6 200 0.02 
0 0.5 0 0 
0 500 1000 0 500 1000 


X [m] X [m] 


Figure 11. Metric components for the comoving modified time metric. The motion 
of the bubble is in the +X direction. The slice is for z = zg, vs = 0.1 c, Amar = 2, 
R = 300 m, o = 0.015 m™t, £o = yo = % = 503 m. 


Figure[l2]ilustrates the energy density, which exhibits similarities to the Alcubierre 
metric. Examining Figure the stress-energy tensor reveals the presence of positive 
pressures in the interior wall of the bubble, in contrast to the negative pressures 
observed in the Alcubierre metric. The modification of the lapse rate results in different 
energy condition violations, as depicted in Figure Lastly, Figure |15| showcases the 
metric scalars, which exhibit similar maximum values but with broader shapes than the 
Alcubierre metric. 
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Figure 12. Energy density for the modified time metric. The motion of the bubble 
is in the +X direction. The slice is for z = zo, vs = 0.1 c, Amaz = 2, R = 300 m, 
o = 0.015 m7, £o = yo = Zo = 503 m. Units are in J/m?. 
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Figure 13. Stress-energy tensor components for the modified time metric. The Tôô 
component is shown in Figure The motion of the bubble is in the +X direction. 
The slice is for z = zo, vs = 0.1 c, Amar = 2, R = 300 m, o = 0.015 m7}, 
£o = Yo = Zo = 503 m. Units are in J/m?. 
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Figure 14. Energy conditions for the modified-time metric. The motion of the bubble 
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Figure 15. Expansion and shear scalars for the modified time metric. The motion 
of the bubble is in the +X direction. The slice is for z = z9, vs = 0.1 c, Amaz = 2, 
R = 300 m, o = 0.015 m™t, £o = yo = % = 503 m. 
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4.4. Lentz-Inspired Metric 


The Lentz metric, introduced in a recent publication [12], presents a dramatic departure 
from the conventional spherically symmetric solutions proposed by Alcubierre and 
others. Unlike its prior warp solutions, the Lentz metric tackles head-on the challenge of 
avoiding violations of the WEC by incorporating shift vectors along multiple directions. 
Lentz’s insight lies in recognizing a limitation of previous warp metrics, where a metric 
featuring a single shift vector component inevitably generates negative energy density 
for Eulerian observers. In contrast, Lentz’s approach involves utilizing additional shift 
components derived from a “potential” described by a linear wave equation. 

A metric inspired by this construction was created in the constant velocity comoving 
frame for analysis in Warp Factory. Instead of the linear wave equation method used 
by Lentz, we construct a piece-wise metric using two shift vector components in the 
X and Y directions. The values and shape of these shift vector sections were designed 
to resemble Figure 1 from [12]. Once these 2D sections were created, the shift vectors 
were smoothed via 2D Gaussian smoothing. Figure [I6] illustrates similar shift vectors 
to those proposed by Lentz. The metric’s parameters are v, = 0.1 c, initial coordinates 
Xp = Yo = 503 m, and a Gaussian smoothing factor of 10 m. 
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Figure 16. Metric components for the comoving Lentz-inspired metric. The motion 
of the bubble is in the +X direction. 


Figure [I 7Jillustrates the evaluation of the energy density wherein the energy density 
remains concentrated primarily at the corners where derivatives of the metric exist. A 
noteworthy distinction in our metric as compared to Lentz’s metric is the distribution 
of Eulerian energy density with ours having regions of both positive and negative 
energy only at the corners of the shift regions. While the energy density predominantly 
approaches zero, violations of the energy conditions emerge consistently at the interfaces 
of the rhomboid structure, as depicted in Figure The cause of this violation lies in 
the non-zero pressure terms at all boundaries of the shift vectors, as seen in Figure 
18| where the energy density is absent. This phenomenon manifests as negative energy 
density to observers in other timelike frames. 
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Figure 17. Energy density for the Lentz-inspired metric. The motion of the bubble 
is in the +X direction. v, = 0.1 c. Units are in J/m?. 


While the reduction of negative energy density in the Eulerian frame is important, 
removing violation requires consideration of the pressure, momenta, and energy density 
relationship, which the Lentz-inspired solution does not manage. From the work here, 
it is not demonstrated that no physical solution exists using this approach, but a key 
challenge has been highlighted for the selection of the gradient boundaries, which must 
generate positive energy density to overcome the non-zero pressure and momentum 
density to make this solution physical. Namely, any physical variant must have positive 
energy along all edges of the rhomboid structure. Lastly, Figure [20] presents the scalars 
for the Lentz-inspired metric. The full construction and evaluation of the Lentz metric 
in Warp Factory was beyond the scope of this paper and will be left for future work. 
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Figure 18. Stress-energy tensor components for the Lentz-inspired metric. The T0 
component is shown in Figure The motion of the bubble is in the +X direction. 
Us = 0.1 c. Units are in J/m?. 
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Figure 19. Energy conditions for the Lentz-inspired metric. The motion of the bubble 
is in the +X direction. vs = 0.1 c and Nobservers = 1000. The minimum value among 
all observers is shown which can be positive. Units are in J/m?. 
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Figure 20. Expansion and shear scalars for the Lentz-inspired metric. The motion of 
the bubble is in the +X direction. vs = 0.1 c. 
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5. Discussion 


5.1. Physicality Results 


For each of the metrics, we performed a full null and timelike vector field analysis of the 
energy conditions, with the results shown in Table|1| Generally, when considering the 
physicality of a warp solution, violations can be attributed to two factors: the presence 
of negative energy density or the occurrence of pressure and momentum density terms 
surpassing the energy density at specific locations in the bubble. Consequently, merely 
ensuring that the energy density is non-negative from the perspective of an Eulerian 
observer is insufficient as a criterion for discussing physicality for a metric. 


Table 1. Comparison of proposed warp drive solution physicality as evaluated using 
Warp Factory. The lapse rate œ means a constant rate whereas & is a modified lapse 
within the drive radius. (X) indicates a violation of the corresponding energy condition 
at at least one point in the space. 


Solution Name aa NEC WEC DEC SEC 
(pseudo-Cartesian) 
Alcubierre [I] =a + BF By x X X X 
bı Õij 
Van Den Broeck [I6] =a Fy Ai X X X X 
pı Vig 
Modifed Time —a? + BP M x X X X 
bı bij 
Lentz-Inspired =a? Bia Bu} x x x x 
Br) Oi 


5.2. Scalars and Momentum Flow 


Expansion is the most frequently discussed scalar in the literature and in the general 
public, particularly the expansion scalar of the Alcubierre metric. However, Natario 
has shown that expansion is not a fundamental feature of warp drives [13]. While not 
a requirement, in the warp drives analyzed in this work, it is observed as a common 
characteristic that there is a contraction in front of the drive and an expansion behind 
it. 

The metrics above all exhibit similar contraction in front and expansion behind 
along the direction of travel, with it occurring along the boundaries of change in the shift 
vector. Likewise, the shear is always positive and similarly exists along the boundaries 
where the shift vector changes. The Lentz-inspired metric is a more complicated case, 
but generally, the interior shift boundaries have a contraction and the exterior has an 
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expansion. The Lentz-inspired metric also has shears along the boundaries of the shift 
vector transitions in a similar manner as the other metrics. 

With Warp Factory, visualizations of the stress-energy tensor reveal more of its 
structure than has been discussed before in the literature. Of particular interest are the 
momentum flux components. For each of the warp drives discussed here, we see a similar 
story of momentum structure, which has a kind of toroidal motion, especially along 
directions that are perpendicular to the axis of motion. This can be seen in the cross- 
section plots of the Alcubierre metric in which the p; stress-energy components change 
sign halfway through the warp bubble wall. This structure can be better visualized 
using the momentum flow analysis which can trace path lines of the momentum flux. 
Animations of the momentum flow lines for Alcubierre, Van Den Broeck, and Modified 
Time are provided as supplemental material to this paper. 


6. Conclusion 


Warp Factory is a specialized toolkit designed for analyzing spacetimes related to warp 
drives. The toolkit addresses the unique analysis needs of the warp community and offers 
comprehensive capabilities for evaluating novel warp metrics’ physicality. Traditionally, 
the analysis of physicality in these metrics has been limited to considering the energy 
density observed by Eulerian observers. However, this approach alone is insufficient 
for a thorough evaluation of the energy conditions. To overcome these limitations and 
enable the exploration of more complex solutions, Warp Factory employs a numerical 
framework. 

By using Warp Factory, researchers can extensively explore a wide range of warp 
metrics, particularly in cases where analytical calculations become impractical. The 
capabilities of Warp Factory are demonstrated in this paper by analyzing a handful of 
popular metrics. Using its numerical approach, Warp Factory provides researchers with 
new insights including 3D visualizations and its ability to do comprehensive analysis of 
stress-energy tensors. A notable finding in this work is a demonstration that metrics like 
the Lentz-inspired solution presented here, which may exhibit solely positive Eulerian 
energy density, can still violate the weak energy condition when examined across all 
timelike observers. 

Using Warp Factory the APL team has also been exploring more warp spacetimes. 
In a separate paper, currently in preparation, a new constant velocity warp drive so- 
lution has been found without energy condition violations [7]. The discovery of this 
exciting result was made possible by utilizing the tools provided by Warp Factory, as 
described in this paper. 


Warp Factory can be found at https: //github.com/NerdsWithAttitudes/WarpFactory 
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Appendix A. Solving for the Tetrad Components 


H . . . . . . 
Assume that ef is a lower triangular matrix (motivated in Section |2.3). 


H j u j j E s 
e= | e& e& e& e| = 
"APT TT (a 
e4 €7 
The metric g,, is given by its components 
gı 92 93 G4 
_ | 92 95 Ge GW 
Juw = 
93 96 gJs 99 
94 97 J9 G10 
Rearranging the tetrad equation: 
gh = heen => ches ow 


Or, in matrix notation, 
E’ gE=n 


0 0 
0 0 
ee 0 
€9 €10 


Solving this system of equations for ef in terms of g,,, gives us: 


-C 
e! = 1 g2 99? + 93 Je J10 + J4 J7 98 — 92.98 910 — 93.97 J9 — 94.96 99 
j C(—D) | g3 97? + 92.96 910 + 94 95 99 — 92 97 99 — 9395 J10 — 94.96 97 
94.96" + 92 97 98 T 93 95 J9 — 92 96 99 — 93 G6 97 — G4 95 98 


ge 99 — 97 98 


el c 
a V BC he 
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(A.3) 


(A.A) 


(A.5) 
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where A, B, C, and D are the leading principle minors of g,, and are given by: 


A= 910 
B= (—go" + 9s 910) 
C = (—95 99° — 98.97" — g10 go” + 2.96 97 99 + 95 98 910) 
D = (92? 99° + 93° gx” + 94° ge” 
— 98 910 92" — 9591093" — 959894 — 9191096 — 91 98 97° — 91.95.99" 
+ 2 910 92 93 Je + 2.98 92 9497 + 295 93 9499 +291 96 97 99 
— 2 92 93 97 99 — 2 92 9496 99 — 2 93 94 96 97 
+91 95 98 910) 


(A.6) 


In order to prevent imaginary values in the transformations, the metric signature, 
and therefore the values of (D,C,B,A), must be (-,+,+,+). 


Appendix B. Error Discussion 
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Figure B1. Plots of errors of Warp Factory’s numerical evaluations for a Schwarzschild 
solution. The error is computed by comparing the numerical T°° value to its true value 
which is zero. The low-resolution grid is 20 m, and the high-resolution grid is 1 m. 
The center of the Schwarzschild metric is at 1000 meters and its radius is 0.01 meters. 
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All 4 types of errors can be seen in Figure The Schwarzschild metric acts 
well to highlight the different types of numerical error of the code since, outside of the 
Schwarzschild radius, the evaluated energy density should be 0. Therefore, any com- 
puted non-zero energy density is due to some type of error. 


Appendix B.1. Edge of Grid Error 


The first, and perhaps most important error from a practical standpoint is edge of grid 
error. The energy tensor is returned in the same grid size as the input metric. At 
the edges of the world, a full finite difference cannot be taken as the points needed to 
evaluate the derivatives at the edges are outside of the world. This can be seen in the 
edges of the plot where all three evaluations have considerable energy density spikes at 
the very bounds of the space. For energy tensor and energy condition plotting, it is 
recommended to exclude a boundary of 2 grid points at the edges of the space. 


Appendix B.2. Finite Difference Discretization Error 


This type of error occurs when the rate of change of the metric components nears 
or exceeds the grid resolution. The error is a direct function of grid resolution. As 
can be seen in the plot, the blue low-resolution line has more error over a wider area 
than the higher-resolution orange line. Even with high resolution, though, the orange 
line does exhibit this type of error as it nears the center of the Schwarzschild metric. 
Metrics with sharp transitions can lead to these finite difference edge effects and may 
give energy/violation that wouldn’t actually exist at those boundaries. 


Appendix B.3. Floating Point Round Off Error 


For double-value precision, there is a maximum dynamic range of numbers that can be 
stored. As can be seen in Figure the errors become very noisy around 107%. The 
double precision limit prohibits storing more discrete values. This precision limit is 
based on the maximum value that exists in the array that is storing the values. 


Appendix B.4. Finite Difference Truncation Error 


This error results from the finite difference not taking a perfect derivative of the function. 
The order of the finite difference (second, fourth, etc.) dictates how accurately the real 
derivative is evaluated. As can be seen in Figure the yellow second-order finite 
difference has higher errors than the orange fourth-order finite difference. 
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Appendix B.5. Alcubierre Numerical Errors 


Figure shows the energy density errors versus the exact analytical expression for 
the energy density of an Alcubierre warp bubble. The radius R, velocity v, and o are 
the same as in Figure |1| The bubble here is centered at 1000 meters and the space is 
extended to show how the errors fall off. 
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Figure B2. Percent error for the numerical evaluation with a 1-meter grid resolution 
of the Alcubierre metric compared to an exact analytical evaluation. The numerical 
error is given as a percent of the maximum energy density of the analytical evaluation. 
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